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Abstract 

We show how the Imphcit Regularization Technique (IRT) can be 
used for the perturbative renormahzation of a simple field theo- 
retical model, generally used as a test theory for new techniques. 
While IRT has been applied successf\illy in many problems involv- 
ing symmetry breaking anomalies and nonabelian gauge groups, 
all at one loop level, this is the first attempt to a generaliza- 
tion of the technique for perturbative renormahzation. We show 
that the overlapping divergent loops can be given a completely 
algebraic treatment. We display the connection between renor- 
mahzation and counterterms in the Lagrangian. The algebraic 
advantages make IRT worth studying for perturbative renormal- 
ization of gauge theories. 

PACS: ll.lOGh, 11.25Db 
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1 Introduction 



Quantum field theoretical predictions of physical quantities should in 
principle be independent of the particular scheme used to renormalize the 
theory. The renormalization program allows to get rid of the singularities 
by redefinition of the parameters in the Lagrangian in a consistent way 
for a renormalizable model. Also, in this process we must make sure 
that the relevant symmetries of the underlying theory are preserved and 
therefore avoid the appearance of spurious anomalies which otherwise 
would have to be controlled order by order in perturbation theory by 
imposing symmetry restoring constraint equations. 

As for the existing regularization schemes whilst for the theories with 
low symmetry content nearly all regulators do a good job, this is not 
the case for most theories of particle interactions in which gauge symme- 
try, supersymmetry (SUSY) and so on play a fundamental role. Dimen- 
sional Regularization (DR) [2j |S! |1I is an efficient and pragmatical method 
which explicitly preserves gauge symmetry. However in the presence of 
dimension specific objects such as 7^ matrices, a suitable generalization 
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of the Dirac algebra must be constructed to be compatible with analytical 
continuation on the space-time dimension. This is the case of the Elec- 
troweak sector of the Standard Model. Since chiral symmetry is broken 
in this case, the corresponding Ward-Slavnov- Taylor identities must be 
imposed order by order, what turns the computations beyond one loop 
order very hard. 

For SUSY theories, the fact that the equality between Bose and Fermi 
degrees of freedom only holds for specific values of the space time dimen- 
sion, SUSY is broken in DR. A naive scheme (Dimensional Reduction) 
in which the field components are left unchanged while the loop inte- 
grals are performed in d dimensions can be shown to be inconsistent, see 
reference. 0. Similar problems arise in Chern-Simons field theories in 
which the Levi-Civitta tensor is the three dimensional analogous of the 
7^ matrix [HllZl- 

A particularly interesting regularization independent framework is 
the Differential Renormalization program pioneered by Freedman et al 
|H]. The basic idea of this scheme is that renormalization comes from 
the fact that products of propagators must be extended to be distribu- 
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tions so that a Fourier transform is well defined. Working in (Euclidean) 
coordinate space one writes the amplitude as a derivative of a distribu- 
tion less divergent at coincident points. The derivatives are understood 
in the sense of distribution theory, i.e. acting formally by parts. The 
amplitudes written in this way are identical to the bare ones for sepa- 
rate points but behave well at coincident points. An intrinsic arbitrary 
scale appears in this process which is used as a Callan-Symanzik renor- 
malization group parameter. The advantage of this method is that it 
works in integer space-time dimension, and it has been shown to yield 
satisfactory results where it was tested [H] [IHl [H] [12] • However no gen- 
eral procedure using Differential Renormalization beyond one loop order, 
such that gauge invariance is automatic, has been constructed yet. 

Recently an essentially regularization independent procedure has been 
advanced (see refs. [ TB] [T7] ) . It presents the same consistency 

as Differential Regularization, working however in momentum space. It 
has been recently shown that the method respects both abelian and non- 
abelian symmetries at a few loops level and the infrared divergence can 
be treated by the same procedure constructed for ultraviolet divergences. 
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without any modification, as has been shown for the gluon self energy 
All these results indicate that the method deserves further investi- 
gation, for example, the question of n-loops and perturbative renormal- 
ization in a theory. Encouraged by the previous results and also aware of 
the difficulties of setting up a consistent scheme for perturbative renor- 
malization where there are already so many of the, so well established, we 
feel for the sake of completeness of the Implicit Regularization Technique 
to take one step in this direction This is precisely the purpose of present 
contribution, Since this is a major task, already pursued by several so- 
phisticated schemes, we should like to start by investigating a simple, 
although non trivial theory (involving, e.g. overlapping divergences) to 
test our method further. If we succeed, it will be most important first 
step to construct a perturbatively n-loop renormalization scheme in the 
spirit and consistency of e.g. BPHZ and others. For our so far very 
humble purposes, we study specifically the (/)g theory. 

The IRT is essentially regularization independent in the sense that a 
specific regulator needs never be used. A convenient identity at level of 
the integrand enables us to rewrite the amplitude as a sum of three types 
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of contributions namely local divergences (basic divergent integrals which 
characterize the divergent structure of theory), nonlocal divergences typ- 
ical of divergent sub-structures contributions) and finite contributions. 
The local divergences obtained in this way are equivalent to those ob- 
tained by performing a Taylor expansion (like in BPHZ) only for prim- 
itive diagrams, with no subdivergences. In the case subdivergences are 
present some more subtle differences between our method and BPHZ 
should be pointed out: a) Taylor expansions modify the original Feyn- 
man integrand; this procedure may therefore violate symmetries. We 
circumvent this problem by means of using a mathematical identity, pre- 
serving thus the original content of the original Feynman amplitude, b) 
In order to classify the subdivergences, BPHZ uses graphic representa- 
tion (forest formula). In our procedure we identify, in terms of integrals, 
divergences which occurred at lower loop orders. This avoids in partic- 
ular the complicated topological graphs structures and substitutes this 
step by an algebraic procedure, (at least in this simple case)from which 
the counterterms appear in a natural systematic way. Specific examples 
are given in the text. Moreover, just like Differential Renormalization 
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arbitrary local terms can be duly parametrized and properly adjusted 
on physical grounds. This is particularly important for finite renormal- 
ization in order to clear the calculation from regularization ambiguities. 
Finally our framework lives in the integer space-time dimension which 
avoids well-known problems with dimension specific theories. 

In order to illustrate our method we study the renormalization of 0g 
theory to n-loop order. 

2 Renormalization by the Implicit Regu- 
larization Technique 

In this section we construct an extension of a technique firstly designed 
for one |13j|l^ and two loop calculations ^SI for performing a n*'*order 
renormalizability proof. 

In order to illustrate the procedure, consider the following divergent 
amplitude, typical of one loop order: 



The symbol A under the integral sign presupposes, as discussed, an im- 
plicit regularization. Now, in order to separate the logarithmic divergence 




— m 
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from the finite part, we use tfie following identity in the factor involving 
the external momentum p: 

1 ^ ^ {-ly {p^ + 2p-ky 

^ (-irv+2p.fcr^ (2) 

(k"^ — m?)^~^^ [{k + pY — m?] 
In the above expression is chosen so that the last term is finite under 
integration over k. Notice also that in the first term in equation the 
external momentum appears only in the numerator and thus after inte- 
gration it can yield at most polynomials in p multiplied by divergences. 
For our present example we need = 0, since we are dealing with a 
logarithmic divergence. We can rewrite (Q) using (j21) as 

/• d'^k 1 r d^k p"^ + 2p ■ k 

" A (27r)4 - m2)2 ~ J (27r)4 [{k + pf - m2](F - m^f ' 

Now only the first of these two integrals is divergent. The others can be 
easily integrated out to yield 

/ = /w(^')- (^^oK,p2) (4) 

where 
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and 



Zo{m'^,p'^) = dz In 



p^z{l — z) — m? 



) 



(6) 



JO 



— 



Note that, since no explicit form for the regulator has been used, one can 
make immediate contact with other regularizations. Details of calcula- 
tions of several one loop amplitudes and their associated Ward identities 
by using this method can be found in [T^ . 

By convenience we divide the diagrams which contribute to a given 
order in two classes: the first which does not contain diagrams which 
possess two point functions as sub divergences and in the second class 
those which do. 

Let us start with the first class of diagrams. To show how the proce- 
dure works it is enough to consider a general Feynman amplitude with one 
external momentum p, one coupling constant A and one mass parameter 
m. We work in the 4-dimensional space-time although the generalization 
to any integer dimension is straightforward. We denote by g a sum of 
internal momenta fcj. The amplitude in question can always be written 
as 
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where 



1 



(8) 



[{p - 



— m- 



and 



1 = number of / structures 



n = number of loops.- 



Note that we have exphcitly separated the terms involving the ex- 
ternal momentum in the denominator, from which nonlocal divergent 
contributions can arise after integration over the internal momenta. The 
structure R{p, q, m, A) contains all other ingredients of the amplitude 
such as coupling constants, results of Dirac traces, and so on. 

For simplicity we adopt the following notation 



r = (ni?)(n/) 



(9) 



where 




(10) 



and 
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As discussed before the source of all possible troubles in the renormaliza- 
tion process will arise from the structure (11/). Our method focus atten- 
tion on these structures. In order to clearly separate finite, local diver- 
gences (whose dependence on the external momenta is only a polynomial) 
from the nonlocal divergences we use a strategy which is completely based 
on the identity 

Define the operator which acts on each structure / in the following 
way 

TO/- ^ I ^P-1^-P' i I 1 (12) 

Qj — m? (g| — m?) \ [{p — qjY — m?] J 

1 I (2p.g,-p^) , {2p.q,~p^f \ 1 \ 

q^j — m? [q^ — m?Y {q^ — m?Y 1 [(P ~ Q'i)^ ~ ""^^l J 

rj.2j ^ 1 ^ {2p.q,-p') ^ {2p.q,-pY 
g| — (g| — m'^y (g| — m^)^ 

, i^p-^-pr i I 1. (14) 

(g| — m^)^ I [{p — q^Y — J 

Note that the action of the operator is equivalent to a Taylor expan- 
sion around zero external momentum where the first terms are kept and 
the rest of the series is resumed, yielding thus a convenient identity. Note 
also that the degree of divergence of the various terms is decreasing. 
The procedure we have in mind consists of applying the operation, in 
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a particular amplitude with superficial degree of divergence to each 
function fj 

r^r = (n/2)n7;^/,.(p,g,,m2)- (is) 

The result of the operation will always have the form 

T^/(p, q, m") = f'%p, q, m') + p'"\p, q, m') ■ (16) 

We define 

f'''{p,q,m') = J2riP,<l,rn')- (17) 

i=0 

Let us exemplify. Take a quadratically divergent amplitude. To each 
contribution of the form 



(P - QjY" - 



\2 ™2 



we associate 



1 

— 777,^ 



f{q,m')^- (18) 

'Yii 'YYi ^ 



(19) 
(20) 



and 



' ' [q^ — m?Y {q^ — rn?Y[{j) — qjY — m?] 

The definitions (|THjl . (fT!I|l . pn|l . pT|l are not unique. It is simply convenient 
for our purposes. Using tliese we rewrite tlie amplitude as a sum of 
various contributions. According to our notation 

T^r = (Hi?) n[/f^(p,g,m2) +/f"(p,g,m2)] ■ (22) 

In this way we can identify three distinct contributions for the amplitude 

T V = ^ fin + ^local "I" ^nonlocal (^3) 

where 

r}.. = (ni?)n//-(p,g,m^)- (24) 

The second contribution contains only local divergences and, for some 
particular {HR) structures, it can contain finite contributions too. It is 
identified as 



r local = (n/2)n//™(p,g, 
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These local divergences correspond to counterterms which are character- 
istic of the order we are renormalizing. For example, they can have the 
form 

The last term in equation ()2H|1 . namely the cross-terms, contain finite 
contributions as well as "nonlocal" divergences. 

^nonlocal = ^ fin "I" ^nonlocal ' (^7) 

These nonlocal divergence contributions will always appear due to the 
divergent subdiagrams (beyond two point functions) contained in the 
graph. As we will show next in a particular example, the renormaliza- 
tion of previous orders will always allow one to cancel these contributions 
if the theory is renormalizable. In the present scheme the result is auto- 
matic and follows from the operation we have just defined, in an algebraic 
manner. There is no need for graphic representations of relevant contri- 
butions, although it is possible. 

The renormalized amplitude say, in n^'^loop order, can therefore be 
defined as 

p(n) _ rpD-p{n) _ j^div{n) _ j^div{n) /nQ\ 
R — -'- local nonlocal \'^°) 
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■pl(n) y2{n) p3(n) 
fin ' fin ' fin 



where the contributions rf™^"-* and ^nontocai contain the counterterms 
typical of order n as well as the counterterms coming from divergent 
subdiagrams of previous order as will become clear in the examples. No- 
tice from the equation above that our framework automatically delivers 
the counterterms 

■pl -pdiv -pdiv ^'OQ^ 

^ CT — local ^ nonlocal l^^J 

and just as in BPHZ, subtracting off the necessary counterterms leaves 
us with the finite part of the amplitude. The main difference between 
our method and BPHZ is that we never modify the original Feynman 
amplitude, since we use an identity at the level of the integrand and 
BPHZ a Taylor expansion. An immediate consequence of this difference 
is that in the present procedure symmetries can always be preserved as 
has been shown in references |13j [n| |15j |16j . 

Now we proceed to evaluate the second class of diagrams, namely 
those which contain two point functions as subdiagrams. Let us call U 
all the two point diagrams contained in a given amplitude P. It is easy 
to see that that they can be factored out inside of the total amplitude in 
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the following sense 

r= n (30) 

all Hj G U 

where TZj stands for the remaining pieces in the amplitude, j characterizes 
a specific two point function, qj is one of the integration momenta (but 
external to E^). Now since the operation T^T is an identity, i.e. T^T — F 
we can define the partially renormalized amplitude (with all two point 
function subdiagrams properly renormalized ) as follows 

f = r + r^^ (31) 

therefore we have 

r^T= n 7^,[5fm^-Afg|] (32) 

all Ej 6 U 

and Tqj. are all counterterms characteristic subdiagrams involving two 
point functions. Sf^rn^ stands for the mass renormalization and for 
the wave function renormalization. Explicit expressions for these objects 
will be given in the following section where a specific example is worked 
out. In order to get the renormalized amplitude of order n from F one 
proceeds in the same way as for diagrams of class one defined above. We 
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thus have 



R local nonlocal 



^in + ^}in + ^%n- (33) 



The whole procedure will become apparent in the concrete example of 
the following section 

3 A(/)g Theory as an example 

Consider the theory Lagrangian, 



^=2 



(dM^)f - ml4{x)] - ^.pUx) (34) 



It is easy to show that a Feynman graph in this theory has the superficial 

degree of divergence D written as 

D^6-2N (35) 

where N is the number of external legs. This means that only Green's 
functions with N < 3 are divergent. For the one-point functions we will 
assume that we can impose the condition (0| |0) = at all orders and 
we will not worry about one-point diagrams. We will just work with 
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the two and three-point Green's functions which possess quadratic and 
logarithmic divergences. 

We will effect the renormalization through the redefinition of the La- 



grangian parameters as: 



(36) 



2 ry 2 



(37) 



-^0 — ZxX 



(3^ 



which allow the Lagrangian to be rewritten as 



£ — Cp + CcT 



(39) 



where 



2 J,2 



A 

3!' 



(40) 



and 



[Z^ - 1) {d,4>f - {Z^Z^ - l)mV'] - {zfZx - 1)^0^ (41) 



At the n*'* order one has 



CcT — >^CT + ^^CT -\- ■ ■ - C,, 



,(n) 
CT 



(42) 



We effect the renormalization at each order imposing the conditions: 
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• Relative to the propagator 



(43) 



and 



p=0 



(44) 



Relative to vertex function 



— iMi^(O) = + finite corrections) 



(45) 



We can rewrite the bare Lagrangian as 



1 r 



(1 + A) {d^(PY - {m' + W)0'J - (1 + B)-<^' . (46) 



A 



in order to identify the renormalization constants 



zl") = 1 + 



(47) 



yin) _ 1 + B^-^ 



(48) 
(49) 



at each order by the imposing renormalization conditions. Since, in prac- 



tice we renormalize each diagram of the given order, the counterterms 
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can be written as 



A('^) = ^a5") (51) 



j 



^Wm^^^^f^m^ (52) 

here a, (b) is the number of three(two) point diagrams which contribute 
to order n . 

At the n*'^ order the inverse propagator function is written as 

Dn'ip') = P^-m^-Ei,^)(p2)- Eg) (/)... 

-5(")m2 + _ s(")(p2) (53) 

and the vertex function as 

-iMR{p,p') = -iX{l + V^'\p,p') + V^'\p,p')... 

+V^''\p,p') + B^^^- (54) 

Using the technique in each diagram contained in the E*^")(p2) and in the 
V^"\p,p') amphtudes we separate the local divergent part and identify 
all divergent substructures. Imposing renormalization conditions we can 
always identify A^"\ S^^'^m'^ and 
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In order to identify the counterterms of the order in question and to 
write the nonlocal ones in terms of divergences of lower orders, showing 
thus that one needs not worry about them, it is convenient to define the 
following functions: 

• Relative to vertex correction counterterms (type j diagrams) 




(n) 



(-^A)^"+^(^)^"7a(m^A) 



local 



(55) 



where 




(56) 



with 




(57) 



For n = 1 



<5 = i, 



(58) 



otherwise 




(59) 
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Notice that what we have defined here are generahzations of the 
simple one-loop logarithmically divergent integral Iiog{m?) which 
we encountered in our one-loop example. 

Relative to all finite contributions to vertex corrections(type j di- 
agram ), which corresponds to the renormalized diagram 



p(n) _ pl(n) p3{n) 
^ R — ^ fin ^ ^ fin 

n p A&t.. 



Relative to the finite contribution, defined in equation ()24|) for the 
overlapping diagrams 



i=l 



In each order there will appear new types of divergent integrals. 
Therefore throughout the text we will define some new divergent in- 
tegrals similar to the ones above (eq. (jKH|) ). These quantities are always 
independent of external momenta. Next we apply the procedure to all 
diagrams up to two loops in order to exemplify how the method works. 
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To n*'* order it suffices to treat four cases, the first related to the vertex 
function and the others to the self-energy, which contain the overlapping 
divergences, two point functions as subdivergences and nested two point 
functions. 

3.1 Three point functions 
3.1.1 The one loop order 

The vertex correction has only one contribution at one loop level whose 
diagram is depicted in figure 1. The corresponding amplitude is 

r = -iv^^\p,p') 

J A (27r)6 - m^) [{p - kf - w?] [{p' - kf - w?] ^ ' 
Using the notation introduced in section 2 we write 

-iV('\p,p')^ / ^^R{k,m',X)f{p,k,m')f{p',k,m') (63) 

with 

According to (IRT) we write, given that the divergence is logarithmic 
and therefore D — 

- iT V(^) (p, p') = + r},„ + (65) 
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(recall that in this case Tj^^ = ^fZiiocai = ) where 



A 



A (27r)6 (A;2 - rn?f 

(66) 



and 



with 



A35(^)(A;,p,p',m2) = R{k,m^,X){ff'''{k,p,m^)f^'''{k,p',m^) 
+ f{k,m')ff^-{k,p',m') 

+ff''^{k,p,m')f{k,m')}- (68) 

Notice that the finite part of this diagram contains the cross-terms 
/o • //in since its integral is finite. 
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3.1.2 The two loop order 

Three diagram types contribute to the vertex correction at two loops. 
The total amplitude can be written as 

- iV'^^^ {p, p') = -3i\//^^ (p, p') - iV^'^^ {p, p') - iV^'^'^ {p, p') ■ (69) 

In this order the counterterms will be identified as 

^(n) ^ 3^(2) ^ ^(2) ^ ^(2) _ ^^Q^ 

The first amplitude — il//^'* (p, p') corresponds to the diagram in figure 
2. This diagram contains a quadratic divergent subdiagram (a first order 
two-point function correction). It can be completely separated in terms 
of the internal momentum ki as mentioned before. Then 

r = -ivi'\p,p') 

r (-a)^(2)^ 

J A (27r)6 {kj - m2)2[(p - A;i)2 - m^][{p' - k^f - m?] 

x{iE«(fcD} (71) 

where iTl^^\k\) is the one loop self-energy amplitude. The one loop 
renormalized self-energy is 

E«(A;?) = Eg)(A;?) + E«(fc?) (72) 
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where 

T.%{kl) = 5(i)m2 - A^^^kl ■ (73) 

Thus the amphtude containing no two point function substructure is 
directly obtained as 



A (27r)6 {kl - m?Y[{j) - A;i)2 - m?][{j)' - k^Y - m?] 
x{zS«(fc?)}. (74) 

All possible nonlocal divergences in this case will be canceled when we 
consider the one loop renormalization. Next we use the IRT for the 
logarithmic divergence. In our notation we obtain 

- zT°f = f f»„, + f + (75) 

with 

/ ^^(^l-^2,m^A) = ^^r^.{^^^S{kl)] = R{k,^m\X). (76) 

The explicit expression for zSSJ^(A;2) will be given in equation ()107|) . Using 
the above expression for -R(/ci, m^, A) 

^tlai = J^-^,Rikum^X)fik,,m')fik,,m') 

= Hxrm!oU^^>^'A) = M"^- (77) 
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We have just defined another logarithmic divergent quantity which is 
characteristic of the two-loop order. Note the explicit appearance of 
coupling constant. This should emphasize the fact that the amplitude 
depends on a two point function subdiagram, which has been properly 
renormalized. All counterterms possessing such type of subdiagram will 
look like this. 

A^ A) ^XgL(j^{>E»> (78) 

The finite part is 

f{h,m')ff^^(p',h,m') 

h, m^)ff'-{p', h, m^)}. (79) 

It is not necessary to give explicit expressions for the finite part and 
therefore we make explicit the divergent contributions only. 

Now we consider the diagram corresponding to the second contribu- 
tion —iV2^\p,p') which belong to class one (figure 3). The amphtude 
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reads 




x/(p, h, m^)fip\ h, m'')f{p', k2, m') (80) 

with 

R{h, h, m\ A) = (^2_^2)(^2_^2)p^_^^)2_^2] • (81) 
Using the IRT we have 

- ir° vi'^ (p, pO = Ttlai + r},, + r„„„,,,,, (82) 

where 



rfjc^a. = (-^A)^(^)«7g,(m^A)=^sf 

/■ d^fci r r/^fc2 {-iXfjif 

J A (2^ A (2^)fi (A;2 - m2)'^(fc2 - m'^f[{ki - - rn?] 

In this type of structure (to all orders) the nonlocal contribution Tnoniocai 



will have the form 

nonlocal fin ~^ nonlocal 



^nonlocal ^ fin ~^ ^nonlocal (^4) 
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and in this case we have 



Note that this term is completely written in terms of one loop contri- 
butions already considered. Therefore it poses no problem to renormal- 
ization. This particular example illustrates a basic difference between 
the present method and others: the subdivergences need not be previ- 
ously identified. They appear algebraically. In cases were it is simple to 
identify the subdivergences , this is not necessarily a great advantage. 
However in higher orders it might become considerably simpler to iden- 
tify all divergent substructures in an algebraic fashion. In fact, as will 
become clear in what follows, the procedure is designed to display all rel- 
evant (to renormalization) subdivergences. The finite contributions can 
be written as 



(85) 




fin 




(86) 



with 
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R{ki, k2,m^, A) 



x{ff'-{h,p,m')ff^-ih,p',m')ff'-ih,p',m') 

+f\h,m')f\h,m')ff^^{h,p',m') 

+f{h,m')ff'-{h,p',m')f^'-{h,p',m') 

+ff^-ih,p,m')fih,m')ff^-{k2,p\m')} 

{-iXf(if{2ki.k2 - kl)E^^\ki,p,p',m^) 
[kl — m?Y[{ki — ^2)^ — m'^] 



The last term in the above equation is obtained by using the operation 
p2|) considering ki as external momentum. This is necessary to identify 
the one loop structure. 

The last two loop diagram —iV^ ip,p') is depicted in the figure 4. 
The corresponding amplitude is 



-iVj^\p,p') = ^— ^{/?(A;i, A;2,m^A)/(p,A;i,m2) 




xf{p',k2,m^)f{p-p',ki - h 



(88) 



with 



R{ki, A;2,m^, A) 



(89) 



{kl-m^){kl 



— m 



mki-k2Y 



Using the IRT we have 



iT'V^'\p,p') = Tt:ai + '^R 



(90) 
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where 

d'^ki r d^k2 1 



Ja (27r)6 7a (27r)6 {kf - m^y{ki - m^y[{ki - k^f - m?Y 
= zA^jg3K,A) = .5f- (91) 

We defined above another logarithmic divergent quantity. This diagram 
type is often called a primitively divergent diagram. Note that there are 
no sub divergences. 
3.1.3 The n-loop order 

As discussed before we now consider only one contribution of each kind. 
The vertex type contribution depicted in figure 5 is the first one. It will 
appear as a substructure of the overlapping self-energy diagram which 
we will also consider. 

The amplitude corresponding to the vertex correction in figure 5 is 

-tvi''\p,p') = {nR){nf) (92) 

with 

where Q is the same function as defined in (jKH|l and The subscript 1 
in Vi"'\p,p') refers to the fact that only one diagram is being considered 
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(type 1 ). The external momentum dependent part (H/) is given 

1 1 



Using the IRT, we get 



r°(n/) = {/°(A;i,m^) + /^-(A;i,p,m2)} 



In the same way we have 

- iT^Vi:"\p,p) = Tiocal + ^]in + ^nonlocal 



Since F^^^ = we can write 



T^div • D (n) 

local — ^-"1 

= (ni?)/°(A;i,m2)n/°(A;„m2) 

= (-^A)^"+^(^)3"/£,(m^A) 

and 



■L nonlocal — fin ~^ ^ nonlocal 
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where 

d' 213"/" 
^nonlocal — (~^'^) '^'^ (0 " 11 L 7nl\6 

X: S^'^H^i, ka,p,p\ m2)T("-«)(A;„+i, .., m")- (99) 

a=l 

Note that the above equation contains subdivergences which have al- 
ready appeared at lower orders and have already been included in the 
Lagrangian. Here we clearly see that the application of the method dis- 
plays all the subdivergencies in an algebraic way. Moreover it stresses 
the inductive character of the method. If we assume that the theory is 
renormalized at {n — 1)*'^ order, the contribution at n*'* order will solely 
depend on structures (finite and divergent) which have already played 
their role at lower orders. Also it is noteworthy that all divergencies and 
finite parts of all previous orders play an important role at n^^ order. 
The above expression clearly displays a difference between this method 
and BPHZ. In BPHZ one would have to deal with the detailed topology 
of the diagram first. Here the counterterms appear algebraically. This 
simplicity may be due to the example we are working with. However, it 
constitutes even in this case a notorious simplification. 
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3.2 Two-point functions 
3.2.1 The one loop order 



The self-energy has only one diagram contribution at one loop level which 
we identify in figure 6. It corresponds to the amplitude 

where 

Using IRT we have 

(r2)zS«(/) = rf-„, + r},„ (102) 

with 



where we have defined 



and 



(104) 
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The finite part is 



^ fin - 



i2 



A2 /■ 



2 J (27r)6 



+ 



/ 



(27r)6(A;2-m2)4[(p-A;)2-m2 



(106) 



The exphcit calculation of the integral in the above equation leads to 



4 (47r)3 

where F{m'^,p'^), for < Arm? ,is given by 



(107) 



2 2N \/4m^ - 



F{m\p') 



and for > Am? , 



2 arctan + tt 



- 2 (108) 



F{m\p^) 



\Jp^ — 4m^ 



vIpI + Vp^ - 4m^y 



-2 ■ (109) 



We now summarize the results obtained so far for one loop the renormal- 



ization, 



and 



(110) 
(111) 



(112) 
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where //Jg\(m^,A) , l!^^^^\m?,K) and /g^)j^(m^,A) are defined in (j^ . 
(fTU^ and (fTMj) . respectively. 
3.2.2 The two loop order 

Two types of diagram contribute to the self energy correction at two 
loops. The total amplitude can be written as 

S(^)(p2) = 2Sf)(/) + sf(p^). (113) 

Here the counterterms to be identified are 

y4(2) = 2Af ) + (114) 

5(2)^2 = 25f)m2 + 5fm2. (115) 

The first amplitude iT!{\p'^) corresponds to the diagram in figure 7. This 
is the same case we have seen in equation 1)711) . Considering the one loop 
renormalization we can write 

Then we apply IRT and obtain 
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in which 



+P'[^€^!^{m\\\l^) - /g,(m^ A^ A)]} (118) 



and 



/a (27r)6 - m2)E 



/,^?(«.^A^A)^/^^^/^^K'(*;)} (120) 



whereas 



'^'V (27r)6 \ (/tf - m2)5 
The second amphtude iE2^''(p2) corresponds to the diagram in figure 



8. It reads 



(122) 



with 



i?(fc„ fc„ m , A) = 2^(^2_^2)(^2_^2)p^_^^)2_^2]} • (123) 
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The same procedure enables us to write 



j.div _ ^ r r(2) / 2 AN 
'-local — cy y-quadxV'^ 1^^) 



^2^2j4|,^M2)(^2^^) - A)]} (124) 



where 



.(2) . 2 A^ _ /■ ^^^1 /■ ^^^^2 

1 

(/ci — m2)2(A;| — m2)2[(A;i — A;2)^ — m 



X -2-1 ri26l 

(fcf-m2)2(fc2-m2)4[(fci-A;2)2-m2] ^ > 



and the finite part coming from this contribution is 



V%^^—{h-r2h-rh) (127) 



with 



' J {2n)^J (27r)6(A;f-m2)3(A;2-m2)3[(A;i-A;2)2-m2] ^ ' 
' J (27r)6y (27r)6(A;?-m2)3(A;2-m2)4pi-A;2)2-m2] ^ ^ 



37 



^nonlocal ^ fin ^nonlocal ' (-'-'^l) 

In terms of functions Q^^\ki,p,'m'^) and T(^)(/cj,m^) we can write 



~ 2 A (27r)6A (27r)6^ 

m^)T*-^''(A;2, rn^) 

+0«(A;2,p,m2)TW(A;i,m2)}- (132) 

Note that T^^^ is (the integrand of a) logarithmic divergence, which, in 
DR would give us 1/e and when multiplied by the remaining pieces of the 
amplitude would produce the celebrated term Inp^/e [TH]. The (other) 
finite contributions are 



and 



= 2 I (2^y (2^® KhM.P.m 



3 {-i^Yiif f d'^h f d'^k2 

~ 1 

2ki.k2 - kl 



2\ 



(133) 



^ ' (27r)6y (27r)6^ 



©'•^-*(A;i,p, m 



2^ 



(^2 — m'^YKki — /c2)^ — rn"^] 
2, 2A;i.A;2 — A;i 



+eW(A;2,p,m- 2/.i./.2-fc^ 



(A;f -m2)3[(A;i - ^2)2 - m2] 
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/^^"(p, h,m^) [/i(fc2, m^) + m^)] 
+ff'^{pM.rn')[h{h,rn') + h{h,m^)]}y (134) 

Summarizing the two loop renormalization constants obtained are, 
= y[^^.X^?(-^A)-/Sl(-^A)] 

+A'[^^?..Cg?^(^^ A^A) -/g2(m^A^A)] (135) 
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<^^'^"^' = -[yC.l(^^A) + A2/SL2(^^A^A)] (136) 



and 



i?^'^ = A1/g,(m^A) + /g3(m^A)] + 3A^Jg2(m^A^A) (137) 



where 



/SiK, A), (m^ A), jg^K, A^ A), cJH"^^ A^ a), /g3(m^ a) 



and 

41l(m^A),/(J,,2("^'•A^A) 

are defined in equations (fnHl) . (f7g|). fT^ . p[|). (fT^ and (fTTn|) 

respectively. 
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3.2.3 The n-loop order 

Let us first consider the overlapping self energy diagram of figure 9. It 
corresponds to the amplitude 



with 



(™) - '""'7''""' {n /, Qc^.*....-^) ■ (139) 

The external momentum dependent part is 
Using the technique we have, as usual 

zT^Ei") = r},„ + Tiocal + r nonlocal (Ml) 

where 

r}.„ = ^'^^^7^^'"" n/ |^e('^)(^„^„...^„;p,m^) (142) 

and 



local — ^ /m "T ^ Zocai 

n 

= (ni?)n{/°(%,rn2) + f (%,m',p) + f (%,rn^p)}•(143) 
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In this case we have F^^^ ^ . The counterterms characteristic of the 
^th QY^QY identified in the equation 

The nonlocal part is 

^nonlocal — ^ fin ^nonlocal i^^^) 

where 



nonlocal 



2 



n-l 



a=l 

n-2 



n-l 
a=l 

a,fe=l J 

From the above equation it becomes clear that the renormalization of the 
self energy to n*^'-order requires all finite functions defined in previous 
self energy diagrams (up to (n — l)*'*-order) as well as all the divergent 
contributions of the three-point functions also to the {n — l)*'*-order. We 
may associate a graphical representation to the equation above and, in 
this way, compare with the BPHZ results. The first term in the equation 
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()145|) contains a sum of — 1 terms comprising a finite functions of the 
type 6 multiplied by the n — a divergent vertex-type functions. The 
second term is the symmetric to the first one (the vertex functions and 
functions swap sides). Finally the last term contains vertex corrections 
to the left and to the right and finite functions in the middle. This can 
be best visualized in the graph which follows (figure 10). 

Notice that in the present procedure no special treatment has been 
given to the overlapping divergencies or to the nested ones, both ap- 
pearing in the self-energy. The reason is that the algebraic procedure 
produces only disjoint divergent contributions. 

In order to complete the renormalization of this theory we will still 
consider two cases, both belonging to the second class defined previously. 
Firstly we consider a specific case where two point functions explicitly 
appear as sub divergences (see figure 11) and the other is an amplitude 
containing an overlapping divergence diagram as substructure (figure 12). 
As we mentioned before, the total integral contains the two point function 
substructures in factorized form. We therefore effect the renormalization 
of the internal propagators directly using the counterterms of order. In 
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this way we immediately obtain f . Let us first consider the case in 
figure 11. This diagram contains subdiagrams involving nested two point 
functions. Following the prescription which display the renormalized 
contributions of previous orders we get 

r = iE^^\p') (146) 

r rl^h 



A (27r)6 {kl - m2)2[(p - A;„)2 - m^] 



with 



d^'K-i 1 



r fl^ k 1 

xi {TlT,^'-\kl) + 5«m2 - A^^^kl) 
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. i-iX?m2l 1^ ■ (149) 

We can substitute the terms in brackets by renormalized function 

^iR'\kl) = (150) 

In this example the two point functions counterterms can be obtained 
for any order n > 1 as 

= nmi)^ i ^.j^. {^^fn-\ki)} (151) 

lAt^ = {-iXf{i)\ (152) 

6^ A (27r)6 {kl - ^ 

A(27r)6(A;2-m2)4^"' ^^-^^ 
Finally we will consider type 3 diagram (which contain a type 1 dia- 
gram) show in figure 12. The corresponding amplitude reads 



(27r)6(A;2 -m2)2[(p-A;„)2-m2 

(153) 

Note that since the structure E^"" ^''(/c^) can be renormalized at (n — 1)*'* 
order (see first example of order n), the n*^ order structure will also be 
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renormalized by a Ffj^^;. Its counterterms have the same form as in the 
previous example. They are 

= (-m^fL§y.J^. (154) 

iA^t^ = i-ixrm (155) 

6^ J A (27r)6 {kl - ^""^^ 
"/a (2nr {kl -m?Y O^Sr 'H^n)}}- 

Note that in this case the three point functions subdiagrams have been 
renormahzed together with the two point subdiagram, since it is con- 
tained in the latter. 

4 Momentum routing independence 

In the examples of the previous sections we have chosen the momentum 
routing in such way as to obtain the simplest form for the final expres- 
sions. Of course, the counterterms so obtained must be independent of 
the particular routing one chooses. In order to exemplify this we consider 
the last example given (type 3 diagram) . One of the possible choices for 
the momentum routing would be to arrange the labels in such a way 
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that external momentum is present in an internal line of the diagram as 
follows 

= (156) 

f d^kn 1 I ■T^"-~^U( —h \^\\ 

Ja (27r)6 {kl - m^)[{p - knf - m^Y 

and also, 

iLt\v^) = I 1^77^ , -UT!-r'\kl)]- 

^ ' ^ 7a (27r)6 (A;2 -m2)2[(p- A;„)2 -m2] I ^ ^ "^J 

(157) 

These two labels must be equivalent, so that the amplitude is momentum 
routing independent as it should. Note that if the amplitude were finite, 
this could immediately be accomplished through a shift p — kn = k'^. 
However, since the amplitude is quadratically divergent, shifts are not 
allowed without the inclusion of surface terms. This point has been 
extensively discussed in our method (see refs. [TTTT ) and a similar 
procedure can cure this problem in the present model. More difficult 
would be theories with gauge symmetries and work along this line is in 
progress. Note that in Dimensional Regularization the problem does not 
appear since shifts are always allowed. 
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5 Conclusion 

We have considered (in the self-energy) all possible complications which 
usually appear in renormalization procedures: overlapping divergences, 
nested divergencies and disjoint ones, all in the same graph at n-loops. 
We have explicitly shown how these problems can be systematically re- 
solved order by order within our technique in a simple example, devoid 
of symmetries abelian or not. However, several generals aspects of the 
method can be learned already from this simple case : 

there will always be a divergent (local) order dependent contribution. 
Also, there will always be a finite contribution composed by the product 
of all finite parts of fj 's. These two structures (divergent and finite) are 
typical of the n^'^ order and poses no problem for renormalization. 

As we have seen in the examples given, the identities we use in the 
integrand leave us then with crossed products of divergent and finite 
contributions. All possible combinations will appear and all of them can 
either be recognized as structures (finite or divergent) already encoun- 
tered in lower order amplitudes or they will give a finite contribution. 

In summary we have presented a completely algebraic method of per- 
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turbative renormalization. All counterterms appear automatically, and 
presents a close correspondence to BPHZ where graphical representations 
are essential. We feel that if our method also works for abelian and non 
abelian theories this could be a major advance. At the one loop level we 
can preserve gauge symmetry if use is made of relations involving diver- 
gent integrals of the same degree of divergence P1],PS1- The difference 
among those integrals is source of both ambiguities and symmetry vio- 
lations .The results we got at one loop are encouraging We are working 
on the application of this method at the two loops level the Quantum 
Electro dynamic (QED). 

We see that the application of this method leads to a relatively simple 
renormalization procedure. There is no need for a graphic representation 
of the relevant contributions. When a diagram has divergent subdia- 
grams the sub divergences need not be previously identified because they 
will appear in an algebraic way. The procedure display all relevant sub- 
divergences. Investigations of advantages of the present method (if at 
all!) over the existing ones remain as our main research interest right 
now. 
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FIGURE CAPTIONS 



Figure 1: One-loop vertex correction —iV^^\p,p') 

Figure 2: The two-loop vertex correction contribution —iVi {p,p') 

Figure 3: The two-loop vertex correction contribution —iV!f^ {p, p') 

Figure 4: The two-loop vertex correction contribution —iV^'^^{p,p') 

Figure 5: The n-loop vertex correction contribution —iVi"'\p^p') 

Figure 6: One-loop self energy iE(^)(p^) 

Figure 7: The two-loop self energy contribution iT!^^\p'^) 

Figure 8: The two-loop self energy contribution iE2^^ (p^) 

Figure 9: The n-loop self energy contribution iT!{\p^) 

Figure 10: Graphic representation of equation (145). 

Figure 11: The n-loop self energy contribution i'E2^\p'^) 

Figure 12: The n-loop self energy contribution i'E^^\p'^) 
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q 




Figure 1: One-loop vertex correction —iV^^\p,p') 




Figure 2: The two-loop vertex correction contribution —iVi {p,p') 




Figure 3: The two-loop vertex correction contribution — {p,p') 
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Figure 5: The n-loop vertex correction contribution —iv}^\p,p') 
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Figure 6: One-loop self energy iS*^^)(p^) 
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Figure 7: The two-loop self energy contribution iHi (p^) 



Figure 8: The two-loop self energy contribution iS2^''(p^) 
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Figure 9: The n-loop self energy contribution iJ]i^\p'^) 
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(n-1) (n-2) (1) 




Figure 10: Graphic representation of equation (145) 
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Figure 11: The n-loop self energy contribution (p^) 
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Figure 12: The n-loop self energy contribution iT,2^\p'^) 
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